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The  United  States  Air  Force  requires  that  Air  Force  aircraft  be  designed  to  be  durable. 
This  requirement  necessitates  an  analytical  demonstration  that  excessive  cracking 
within  the  airframe  will  not  occur  during  the  aircraft's  design  service  life.  In  order 
to  predict  the  time  at  which  excessive  cracking  occurs,  an  analysis  is  needed  which  is 
capable  of  predicting  the  distribution  of  crack  sizes  within  the  airframe  at  any  point 
in  time.  Such  an  analysis  was  recently  developed  and  is  presented  in  this  paper.  The 
durability  analysis  is  based  on  a fracture  mechanics  philosophy,  combining  a probabilis- 
tic format  with  a deterministic  crack  growth  rate  relationship.  Essential  elements  of 
the  methodology  are  presented,  with  emphasis  on  the  statistical  representation  of  the 
initial  fatigue  quality  of  the  structure.  The  accuracy  of  the  durability  analysis  is 
demonstrated  by  correlating  analytical  predictions  with  experimental  results  of 
a fighter  full-scale  test  article  as  well  as  complex-splice  specimens  subjected  to  a 
bomber  load  spectrum. 


I INTRODUCTION 


The  United  States  Air  Force  structural  integrity  requirements  for  metallic  air- 
frames are  specified  m Reference  1.  Of  particular  importance  are  the  damage  tolerance 
(Refs.  2 and  3)  and  durability  (Refs.  3 and  4)  requirements.  The  '•’amage  tolerance 
requirements  ensure  aircraft  safety  while  the  durability  requireme.  :s  minimize  structural 
maintenance  costs  and  functional  impairment  problems.  Both  sets  of  requirements  are 
necessary  to  ensure  the  operational  readiness  of  Air  Force  aircraft. 

This  paper  addresses  the  durability  aspects  of  structural  integrity.  The  durability 
damage  mode  considered  is  fatigue  cracking  in  fastener  holes.  This  was  found  to  be  a 
very  prevalent  form  of  degradation  in  aircraft  structures  (Rat.  5).  Aircraft  structural 
durability  involves  many  fastener  holes  in  various  components  which  arc  susceptible  to 
cracking  in  service.  The  associated  structural  maintenance  costs  arc  proportional  to 
the  number  of  fastener  holes  requiring  repair.  Therefore,  to  assess  the  durability  of 
the  structure  or  the  extent  of  damage  as  a fmction  of  time,  the  entire  population  of 
fastener  holes  must  be  considered.  Thus,  a statistical  approach  is  best  suited  for 
quantifying  the  extent  of  damage  as  a function  of  time. 

Various  aspects  of  structural  durability  have  been  considered  m the  literature 
(Refs.  o-B.7).  Structural  durability  can  be  defined  in  many  different  ways.  The  most 
appropriate  definition  is  dependent  cn  the  particular  aircraft  considered.  The  crack 
sizes  oi  interest  are  a function  of  the  definition  used.  This  paper  considers  relatively 
small  subcntical  crack  sizes  which  often  affect  durability.  For  example,  0.76  mm  - 
1.27  mm  (C.03  inch  - C.05  inch)  radial  cracks  m fastener  holes  can  affect  functional 
impairment,  structural  maintenance  requirements,  and  life-cycle  costs.  Such  cracks  do 
not  pose  an  immediate  safety  problem.  However,  if  the  fastener  holes  containing  such 
cracks  are  not  repaired,  economical  repairs  cannot  be  made  when  these  cracks  exceed  a 
limiting  crack  size.  For  example,  a 0.76  mm  - 1.27  nun  radial  crack  in  a fastener  hole 
can  be  cleaned  up  by  reaming  the  hole  to  the  next  nominal  hole  size.  When  the  crack 
sizes  exceecs  this  economical  repair  limit,  excessive  maintenance  and  repair  costs  can 
occur.  Hence,  an  analysis  is  needed  for  predicting  the  extent  of  damage  present  at  any 
particular  point  in  time.  Such  an  analysis  was  recently  developed  (Refj.  11-16  and 
18-27)  and  evaluated  for  load  transfer  coupon  specimens  (Refs.  11,  12,  15,  21,  22,  and  27). 
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/0~  This  paper  presents  the  recently  developed  durability  analysis  methodology.  The 

methodology  is  based  on  a probabilistic  fracture  mechanics  approach.  The  initial 
fatigue  quality  of  a structure  is  represented  in  a statistical  manner  by  a distribution 
of  equivalent  initial  flaw  sizes.  These  equivalent  initial  flaw  sizes  are  grown 
forward  in  time  using  a deterministic  crack  growth  rate  relationship.  An  evaluation  is 
made  of  the  accuracy  of  the  analysis  by  correlating  analytical  predictions  with  test 
data  for  a fighter  full-scale  test  article  (Refs.  21  and  22)  and  complex  splice  speci- 
mens subjected  to  a bomber  load  spectrum  (Ref.  27)  . 

II  DURABILITY  DESIGN  REQUIREMENTS 

Air  Force  durability  design  requirements  for  metallic  airframes  are  presented  in 
References  1,  3 and  4.  According  to  these  requirements,  the  airframe  ruu^t  be  designed 
to  have  an  economic  life  greater  than  the  design  service  life.  Furthermore,  the 
economic  life  must  be  demonstrated  by  analysis  and  test.-  The  economic  life  analysis 
must  account  for  the  effects  of  initial  quality  variations,  material  property  vari- 
ations and  the  design  loads/environments. 

The  economic  life  of  a structure  is  currently  defined  in  only  qualitative  terms: 
"...  tne  occurrence  of  widespread  damage  which  is  uneconomical  to  repair  and,  if  not 
repaired,  could  cause  functional  problems  affecting  operational  readiness"  (Ref.  1). 
There  is  no  universal  quantitative  definition  of  "widespread  damage"  or  acceptable 
structural  maintenance  cost  limits.  Such  limits  are  currently  determined  by  the 
contractor  and  the  Air  Force  for  the  particular  aircraft  of  interest.  Durability 
compliance  standards  are  defined  based  on  the  results  of  the  full-scale  durability  test 
article. 

III  DURABILITY  ANALYSIS  CRITERIA 


i 


A.  Durability  Critical  Parts  Criteria 

Criteria  must  be  developed  for  determining  which  parts  of  an  aircraft  are 
durability  critical  (i.e.,  which  parts  must  be  designed  to  meet  the  durability  design 
requirements).  The  durability  critical  parts  criteria  vary  from  aircraft  to  aircraft. 
They  are  especially  dependent  on  the  definition  of  economic  life  for  the  particular 
aircraft  involved.  A more  detailed  discussion  of  the  durability  critical  parts  criteria 
is  presented  elsewhere  (Ref.  12).  A typical  flow  diagram  for  selecting  which  parts  are 
durability  critical  is  presented  m F i.  1.  In  Fig.  1,  durability  refers  to  the 
ability  of  an  airframe  to  resist  crac)  ng  whereas  damage  tolerance  refers  to  the 
ability  of  an  airframe  to  resist  failure  due  to  the  presence  of  such  cracks. 
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Figure  1 - Flow  diagram  for  selecting  durability  critical  parts 


B.  Economic  Life  Criteria 


Criteria  must  be  developed  for  determining  the  economic  life  of  the  particular 
aircraft  of  interest.  Similar  to  the  durability  critical  parts  criteria,  economic  life 
criteria  vary  from  aircraft  to  aircraft.  They  may  be  based  on  fastener  hole  repair 
(e.g.,  reaming  the  damaged  fastener  hole  to  the  next  nominal  hole  size),  functional 
impairment  (e.g.,  fuel  leakage),  residual  strength,  etc.  Two  promising  analytical 
formats  for  quantifying  the  economic  life  of  an  airframe  are  (1)  the  probability  of 
crack  exceedance,  and  (2)  cost  ratio:  repair  cost/replacement  cost.  Both  formats 

require  a durability  analysis  methodology  capable  of  quantifying  the  extent  of  aircraft 
structural  damage  as  a function  of  service  time.  For  example,  assume  the  economic  life 
criteria  are  based  on  the  number  of  fastener  holes  which  cannot  be  economically  repaired 
(i.e.,  number  of  fastener  holes  with  crack  sizes  equal  to  or  greater  than  specified 
size  X^) . Then  an  analytical  format  for  quantifying  economic  life  is  presented  in  Fig. 
2.  InFig.  2,  P is  the  exceedance  probability.  More  detailed  discussions  of  economic 
life  criteria  are  presented  elsewhere  (Refs.  11,  12,  15,  and  18). 


Figure  2 - Analytical  format  for  economic  life 


IV  DURABILITY  ANALYSIS  METHODOLOGY 
A.  General  Description 

The  basic  objective  of  the  durability  analysis  methodology  is  to  quantify  the 
extent  of  damage  as  a function  of  service  time  for  a given  aircraft.  The  extent  of 
damage  is  measured  by  the  number  of  structural  details  (e.g.,  fastener  holes,  cut-outs, 
fillets,  lugs,  etc.)  expected  to  have  a crack  size  greater  than  a specified  size  at  a 
given  service  time.  Hence,  the  extent  of  damage  is  represented  by  the  probability  of 
crack  exceedance.  The  durabilJty  analysis  results  provide  a quantitative  description 
of  the  extent  of  damage  and  a basis  for  analytically  assuring  that  the  economic  life  of 
the  structure  will  exceed  the  design  service  life. 

The  durability  analysis  includes  two  essential  steps:  (1)  the  quantification 
of  the  initial  fatigue  quality  of  the  structural  details  considered,  and  (2)  the 
prediction  of  the  probability  of  crack  exceedance  based  on  the  initial  fatigue  quality 
and  tne  applicable  design  conditions  (e.g.,  load  spectrum,  stress  levels,  percent  load 
transfer,  etc. ) . 


B.  Initial  Fatigue  Quality 

1.  Initial  Fatigue  Quality  Description 

Initial  fatigue  quality  (IFQ)  defines  the  initial  manufactured  state  of 
a structural  detail  or  details  with  respect  to  initial  flaws  in  a part,  component  or 
airframe  prior  to  service.  The  IFQ  for  a group  of  replicate  details  is  represented  by 
an  equivalent  initial  flaw  size  (EIFS)  distribution.  An  equivalent  initial  flaw  is  a 
hypothetical  crack  assumed  to  exist  in  a detail  prior  to  service.  An  equivalent 
initial  flaw  size  is  the  initial  size  of  a hypothetical  crack  which  would  result  in  an 
actual  crack  size  at  an  actual  point  in  time.  An  arbitrary  crack  size,  a.,  is  selected 
which  can  be  readily  detected  or  which  can  be  reliably  observed  fractographically 
following  testing.  The  time  required  for  an  initial  defect,  of  whatever  type,  to 
become  a fatigue  crack  of  size  aQ  is  defined  as  the  time-to-crack-mitiation  (TTCI) . 
Test  results  of  TTCI  and  crack  growth  rates  using  coupon  specimens  are  employed  to 
define  the  EIFS  distribution.  A conceptual  description  of  the  initial  fatigue  quality 
model  is  shown  in  Fig.  3 (Refs.  11,  12,  15,  21).  Once  the  EIFS  distribution  has  been 
defined  for  a group  of  details,  a deterministic  crack  growth  analysis  is  used  to  grow 
the  entire  EIFS  population  to  any  service  time  T;thus  determining  the  time-varying  crack 
size  distribution  in  service  (Fig.  4). 
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Figure  3 - Initial  fatigue  quality  model 


2.  Initial  Fatigue  Quality  Model 

A prototype  initial  fatigue  quality  (IFQ)  model  has  been  previously 
described  for  quantifying  the  EIFS  distribution  for  structural  details,  such  as  fastener 
holes  (Refs.  11,  15,  and  16).  Such  a model  is  conceptually  described  in  Fig.  3.  IFQ 
model  refinements  and  EIFS  distributions  for  two  model  variations  are  summarized  in  the 
following. 

The  TTCI  distribution  for  coupon  specimens  is  represented  by  the  three-parameter 
Weibull  distribution,  F„(t)  (Refs.  11,  15  and  16). 


FT(t)  = P [T  < t]  = 1 - exp  [-( 


-1  ' 1 , 


where  T = TTCI , 


= shape  parameter, 


The  three  Weibull  parameters  (i.e.,  , £.,  ) 

coupon  specimens  or  other  suitable  test  results 


,'0  = scale  parameter  and  e = lower  bound  TTCI. 
0,  £q,  ) are  determined  from  fractography  for 
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Figure  4 - Growth  of  EIFS  distribution  as  function  of  time 


The  crack  growth  rate  in  the  small  crack  region  is  assumed  to  be  of  the  following 

form 

Hr1  =Q0[a(t)]b  (2) 

where  QQ  and  b are  parameters  depending  on  loading  spectra,  structural  and  material 
properties,  etc.;  a(t)  is  the  crack  size  at  time  t.  Other  functional  forms  for  the 
cracK  growth  rate  could  also  be  used.  Equation  2 is  used  because  of  its  simplicity  and 
general  applicability  for  matching  the  crack  growth  data.  A crack  growth  rate  equation, 
such  as  Eq.  2,  is  used  to  obtain  the  EIFS  distribution  through  a transformation  of  the 
TTCI  distribution;  hence  both  EIFS  and  TTCI  are  statistically  compatible.  EIFS  distri- 
butions are  obtained  using  Eq.  2 for  b 4-  1 (case  I)  and  b = 1 (case  II).  The  resulting 
EIFS  distribution  for  both  cases  are  described  below. 

a.  Case  I (b  # 1) 


im 


kViN 


'V-\ 


Integrating  Eq.  2 from  t = 0 to  t = T,  the  relationship  between  the 
initial  crack  size,  a(0),  and  the  reference  crack  size,  aA  = a(T)  for  TTCI  (flicht 
hours) , is  obtained,  0 

EIFS  = a(0)  = (a"c  + cQ0T]~1/c  (3) 

where  c = b - I. 


The  EIFS  cumulative  distribution,  F ,A. 
as  follows:  a(0) 


Fa(0) (x)  = exp 


x“c  - a.c  - cQ^e"1  ° 


-]  } 


(x)  , is  obtained  by  combining  Eqs.  1 and  3 


x < x 
— u 


,(4) 


u 


1; 


x > x 
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JO  - (*  where 


x 

u 


[a: 


+ cQ0e]"1/C 


(5) 


is  the  upper  bound  of  the  EIFS  distribution.  The  lower  bound  of  the  TTCI  distribution, 
e,  for  a given  reference  size,  a^,  is  obtained  from  Eq.  9 as  follows: 


e=-db[x  uC'a0C];  a0  — xu 


(6) 


Equation  4 for  F ....  (x)  can  be  simplified  by  substituting  the  expression  for  CQQe  from 
Eq.  5 into  Eq.  4aas  follows: 


Fa (0) (x)  = eXp 


-c  -c  “ 

- xu 

*0  J 

cQoeo 

> 

) 

= 1; 

b.  Case  II  (b  = 1) 


x < x 

— u 


X > X 

- u 


(7) 


Integrating  Eq.  2 from  t - 0 to  t = T and  considering  b = J,  one  obtains  the 
relationship  between  the  initial  c.ack  size,  a(0),  and  the  reference  crack  size,  a.,  as 
follows  v 


EIFS  = a(0)  = <r.  exp  (-QQT) 


(8) 


Tne  ritrs  cumulative  distribution,  F ...  (x)  , is  obtained  by  combining  Kqs.  1 ana  d 

a (u) 


as  follows: 


Fa(o),<x)  = exp  i- 


fitn  (x  /x)"!  °) 


Q0S0  J 


C < x < x 


(9) 


- 1; 


x > x 


where  is  the  upper  bound  of  the  EIFS  distribution: 


x = 
u 


aQ  exp  (~Q0£) 


(10) 


The  lower  bound  of  TTCI  distribution,  c,  for  a given  reference  crack  size,  a.,  is 
obtained  from  Eq.  10  as  follows: 


in  (aQ/xu)-  aQ  > xu 


(11  > 


When,  each  test  specimen  consists  of  Z fastener  holes  equally  stressed  and  trac tography 
results  are  taken  only  for  the  fastener  hole  m each  specimen  which  has  the  largest 
crack  siz1-*,  the  exponential  exponents  of  Eqs.  1,  4,  7 and  9 should  be  multiplied  by  a 
factor  or  1/5  (Ref.  22). 


Generic  EIFS  and  Discussion  of  IFQ  Model 


Intuitively,  the  EIFS  cumulative  distribution,  F ...  (x) , should  be  only  a func- 
tion of  the.  matex'ia!  and  the  manufacturing/fabrication  processes . As  such  F ...  (x) 

(Eqs.  4,  7 and  9)  should  be  mdepen  lent  of  loading  spectra,  stress  level,  percent  shear 
load  transfer  through  the  fasteners,  etc.  If  the  same  EIFS  cumulative  distribution  is 
valid  for  replicate  fastener  holes  under  different  design  conditions  (e.g..  loading 
spectra,  stress  level,  i load  transfer,  etc.),  tnen  the  resulting  F ...  (x)  is  said  ':o 

be  "generic".  If  r (x)  is  generic,  then  th.  crack  growth  damage  accumulation  can 

be  calculates  analytically  or  numerically  for  different  design  conditions  using  the 
EIFS  distribution  (Eqs.  4,  7 oi  9). 


The  TFQ  model  parameters,  iQ,  3^,  :,  Qq,  and  t depend  on  the  f ractographic  results. 
Therefore,  these  parameters  depend  on  the  conditions  used  to  generate  the  crack  initia- 
tion and  crack  growth  data  (e.g.,  matarial,  loading  spectra,  stress  level,  etc).  A 
basic  premise  of  the  durability  an  .lysxS  method  proposed  is  that  F ...  \x)  can  be  det- 

ermined from  a given  fiactography  lata  -.et  or  sets  economically  ana  tne  resulting  EIFS 
distribution  can  be  used  for  era<_k  xc 'edance  predictions  for  different  load-ng  spectra, 
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stress  levels,  % load  transfer,  etc.  Encouraging  results  have  been  obtained  to  date 
which  suggest  that  such  a basic  precise  appears  to  be  promising.  However,  further 
study  is  required  to  evaluate  the  IFQ  distributions  using  available  f ractographic 
data  generated  in  Ref.  23  and  to  assess  the  accuracy  of  the  crack  exceedance  pred lo- 
tions, p(i,x),  under  different  design  conditions. 

For  simplicity,  suppose  two  sets  of  replicate  specimens  are  tested  using  the  same 
loading  spectrum  but  different  stress  levels.  Using  the  f ractographic  results  for  each 
data  set,  the  respective  T'i'CIs  for  a given  a.  can  be  determined  for  each  data  set  as 
illustrated  m Fig.  5.  For  the  EIFS  cumulative  distribution,  F ...  (x) , to  be  "generic" 

the  TTCI  distributions,  F_(t),  for  data  sets  1 and  2 should  transform  into  the  same 
EIFS  distribution.  From  Eq.  9,  the  necessary  conditions  for  a generic  EIFS  cumuiative 
distribution  are  for  «q  and  Qq8q  to  be  constants. 


Id 


Figure  5 - Generic  EIFS  condition 


& 


Investigations  have  snown  that  is  a material  constant  for  a given  type  of 
fastener  hole  and  the  product  Q.Bq  appears  to  be  the  common  denominator  for  linking 
different  f ractographic  data  setsutogether  on  a common  baseline.  Although  results  to 
date  are  very  encouraging,  futher  research  is  required  to  evaluate  and  compare  Q.6q 
values  for  different  f ractographic  data  sets. 

The  EIFS  distnbucio  : for  case  I (b  # 1)  and  case  II  (b  = 1)  are  giver,  by  Eqs.  7 
and  9,  respectively.  In  both  cases,  the  EIFS  cumulative  distribution,  F ... (x) , is 
independent  of  the  reference  crack  size,  a.,  used  to  define  the  TfCI  values.  This  is 
an  important  attribute  of  the  IFQ  distribution.  The  upper  bound  EIFS  value  in  the  IFQ 
distribution,  X , is  either  selected  by  the  user  or  computed  from  Eq.  5 or  10.  The 
other  parameter^  (i.e.,  a.,  Bq,  e,  C.  and  b)  are  determined  from  f ractographic  results 
for  coupon  specimens  or  other  suitable  test  results. 

In  previous  work,  Case  I (b  ^ 1)  has  been  considered  (Refs.  II,  12,  16,  21,  and 
22).  Fractographic  results  for  protuding  head  (Ref.  33)  and  countersunk  head  (Ref.  23) 
fasteners  resulted  in  b values  <1.  When  b isi.  1.0,  it  is  poss.  ble  to  obtain  EIFS 
values  £0  using  the  original  IFQ  model.  Mathematically,  however,  the  original  IFQ 
model  can  handle  both  positive  and  negative  EIFS  values,  since  the  IFQ  model  is  simply 
a "mathematical  tool"  for  predicting  the  probability  of  crack  ejxceedance.  An  explanation 
of  the  negative  EIFS  issue  m terms  of  the  original  IFQ  model  and  the  crack  initiation 
process  is  given  m Ref.  22. 


'.■■here  are  certain  advantages  and  disadvantages  associated  with  the  use  of  Case  I 
(b  ^ I)  and  Case  II  (b  = 1)  in  Fq.  2.  For  example,  the  EIFS  master  curve  for  Case  I 
v.’ith  two  parameters  Qn  and  b j5  1 generally  fit  the  selected  fractographic  crack  size 
range  better  than  Case  II  with  a single  parameter  (Q,.,  b = 1).  However,  the  resulting 
Qq  and  b values  must  be  on  a comparable  baseline  when  different  fractography  data  sets 
are  considered.  Since  the  resulting  Q.  aid  b values  affect  the  IFQ  distribution,  they 
must  be  consistently  defined.  In  Eq.  2,  Q.  and  b are  shown  to  be  strongly  correlated 
parameters  (Ref.  19) . Hence,  for  a given  B there  is  likely  to  be  a corresponding  Q„ 
and  vice  versa.  Thus,  for  consistent  IFQ  results,  the  same  b value  may  be  determined 
using  pooled  fractography  results  for  different  data  sets. 


When  the  one  parameter  form  of  Case  II  is  used  in  Eq.  2 (i.e. . Q_,  b = 

"Q  " value  for  one  fractography  data  set  is  already  comparable  with  tne  "Oq 
another  data  set.  Nevertheless,  whichever  form  of  Eq.  2 is  used  (i.e.,  Q„;  b 
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b = 1)  the  resulting  F 


(x) 


distribution.  Consequently,  as  long  as  the  resulting  F 


lcnever  xorm  oi  r.q.  t is  used  (i.e.,  y. , 
will  be  statistically  compatible  with  th« 


manner,  the  same  crack  exceedance  prediction  will  be  o'. 
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1) , the 
value  for 
* 1 or 
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(x)  is  used  in  a consistent 
eJ. 


C.  Durability  Analysis  "rocedures 


The  durability  analysis  procedures,  described  and  discussed  in  detail  elsewhere 
(Refs.  15,  21  and  22),  are  summarized  below  for  Case  II  (b  - 1).  Similar  expressions 
can  be  developed  for  Case  I (b  = 1)  using  the  same  procedures. 

(1)  Divide  the  durability  component  into  m stress  regions  where  the  maximum 
stress  m each  region  may  be  reasonably  assumed  to  be  equal  for  every  location  or 
detail  (e.g.,  fastener  hole). 


(2)  Use  the  model  shown  in  Fig.  3 and  suitable  fractography  results  to  define  the 

EIFS  distribution  expressed  in  Eq.  9 ript-rrminp  n and  Q « (Ref,  22).  The  x selected 
should  be  consistent  with  Eq.  10.  u 

(3)  Determine  for  each  stress  region,  the  corresponding  EIFS  value,  y^.  (tj,  which 
grows  to  crack  size  x.  at  service  time  t as  illustrated  in  Fig.  4 (Ref.  22)7  If 
applicable  fractography  data  are  available  for  different  stress  levels  and  fractography 
data  pooling  procedures  are  used,  the  crack  growth  rate  expression  in  Eq.  12,  where  b^  = 1 
can  be  integrated  from  a(0)  -=  y^fx)  to  a(x)  - x^  to  obtain  y^tx)  in  Eq.  13. 


ir1  = QiCa(t)j  1 


(12) 


yli  (t)  = X1  exp  (-QiT) 


(13) 


in  which 

Q.  = £oY 

i 

is  assumed  to  be  a power  function  of  the  maximum  applied  stress  a, 
constants  ti  be  determined  from  available  fractography  data. 


(14) 


and  £ and  y are 


If  applicable  fractographic  results  are  not  available  for  the  desired  design 
conditic..s  (e.g.,  load  spectra,  % load  transfer,  stress  level,  etc.),  an  analytical 
crack  growth  program  (e.g..  Ref.  36)  can  be  used  to  generate  a "service  crack  growth 
master  curve"  to  determine  y.  (x)  for  a given  x,  (Refs.  11  and  15) . Wr.en  an  analytical 
crack  growth  program  is  used,  y,  (x)  must  be  determined  in  a manner  which  is  consistent 
witn  that  used  to  determine  theAEIFS  distrib jtion  from  the  fractographic  test  data. 

(4)  Compute  the  probability  of  crack  exceedance  for  each  stress  region,  i.e., 
p(i,x)  = P [a  (x ) > xx)  = 1 - Fa  (Q)  fyxi  (t)  ] , using  Eq.  9. 


p(i,x)  = 1 - exp 


1 

en(xu/yli(T) 

a 0 

i 

1 

o 

o 
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1 

0 < y.  ■ < x 
Jli  - u 


(15) 


P (if  t ) 
in  which  y 


0; 


yli(A> 


> X 

- u 


based  on  tii£  number  of  fastener  holes  per  specimen  used  in  the  fractography  data  base. 


(t)  is  given  by  Eq.  13,  and  i ir  the  scaling  fact 


■r.a  orcvicusiv 


(5)  The  average  number  of  details  N(i,  i) , and  the  standard  deviation  o(i,  x)  m 
tiie  ith  stress  region  with  a crack  size  greater  than  x.  at  service  time  x are  determined 
using  the  binomial  distribution  and  are  expressed  as  follows: 

N(i,x)  = NiP(i,x)  <1S* 


= {N,  p ( i , x ) ( 1 - p(i,x)]}1,/2 


(17) 


in  which  N . denotes  the  total  number  of  details  in  the  ith  stress  region.  The  average 
number  of__aetails  with  a crack  si2e  exceeding  x^  at  the  service  time  t for  m stress 
regions,  L (i),  and  its  standard  deviation,  a T(x),  can  be  computed  using  Eqs.  18  and 
19.  L 


m 

L(T)  = X!  N (i,  T ) 
i=l 


°l(t) 


(18) 


(19) 


Equations  18  and  19  can  be  used  to  quantify  the  extent  of  damage  for  a single  detail,  a 
group  of  details,  a part,  a component,  or  an  airframe.  Upper  and  lower  bounds  for  the 
prediction  can  be  estimated  using  L{z)  ±Z c^fx),  where  Z is  the  number  of  standard 
deviations,  J^tx),  from  the  mean,  L(x).  Equations  16  through  19  are  valid  if  cracks  in 
each  detail  are  relatively  small  and  the  growth  of  the  largest  crack  in  each  detail  is 
not  affected  by  cracks  in  neighboring  details.  Hence,  the  crack  growth  accumulation 
for  each  detail  is  statistically  independent  (Refs.  11,  12,  15  and  21). 

V.  DURABILITY  ANALYSIS  DEMONSTRATIONS 


A durability  analysis  of  the  lower  wing  skins  of  afighteris  presented.  Dura- 
bility analysis  results  for  complex-splice  specimens  subjected  to  a bomber  load  spectrum 
are  also  presented.  Both  analyses  are  correlated  with  test  data. 

A.  Fighter  Lower  Wing  Skins 

A durability  analysis  of  the  lower  wing  skins  of  a fighter  durability 
test  article  is  presented  to  illustrate  the  methodology  described.  Analytical  pre- 
dictions of  the  extent  of  damage  in  each  wing  skin  are  presented  in  various  formats, 
and  results  are  compared  with  observations  from  the  tear-down  inspection  of  the  fighter 
durability  test  article. 

The  fighter  durability  test  article  was  tested  to  16,000  flight  hours  (equivalent 
to  2 service  lives)  using  a 500-hour  block  spectrum.  Each  wing  received  the  same 
loading.  Following  the  test,  all  fastener  holes  in  the  lower  wing  skins  were  inspected 
using  eddy  current  techniques.  Fastener  holes  with  crack  indications  were  confirmed  by 
f ractographic  evaluation.  The  right  hand  and  left  hand  lower  wing  skins  were  ^aund  to 
have  twenty  six  and  seven  fastener  holes,  respectively,  with  a crack  size  ^0.76  mm 
(0.03  inch)  at  t = 16,000  flight  hours. 

A preliminary  durability  analysis  for  che  fighter  lower  wing  skins  was  presented 
in  Ref.  21.  The  preliminary  analysis  reflected:  (1)  fastener  hole  IFQ  based  on 

fractographic  results  for  protruding  head  fasteners,  (2)  crack  growth  rates  for  the  IFQ 
model  based  on  Eq.  2 (b  # 1) , (3)  three-parameter  Weibull  distribution  used  in  the  IFQ 
model,  (4)  model  parameters  based  on  a single  data  set  (one  stress  level,  400-hour 
block  spectrum,  Ref.  33),  (5)  three  stress  regions  considered  for  the  lower  wing  skin, 
and  (6)  an  analytical  crack  growth  program  (Ref.  36)  and  the  500-hour  block  spectrum 
were  used  to  define  the  "service  crack  growth  master  curve"  for  each  stress  regirn. 

Essential  features  of  the  present  analysis  are:  (1)  fractographic  results  for 

countersunk  fasteners  used  to  quantify  IFQ  (countersunk  fasteners  were  used  on  the 
fighter  durability  test  article) , (2)  crack  growth  rates  for  the  IFQ  model  based  on  Eq. 

2 (b  = 1) , (3)  three-parameter  Weibull  distribution  used  in  the  IFQ  model,  (4)  mcdel 
parameters  based  on  three  different  data  sets  (three  stress  levels,  400-hour  block 
spectrum)  (5)  lower  wing  skin  divided  into  10  stress  regions,  and  (6)  crack  growth  rate 
parameter  Q defined  for  each  stress  region  as  a function  of  stress  level  a,  determined 
from  available  fractography  data.  There  were  no  significant  differences  in  the  400-hour 
and  500-hour  spectra. 

The  fighter  lower  wing  skin  was  divided  into  ten  stress  regions  as  shown  in  Fig. 

6.  Applicable  stress  levels  and  the  corresponding  number  or  fastener  holes  in  each 
stress  region  are  shown  in  T.  ble  1.  The  stress  levels  for  Zones  I-IV  were  determined 
using  strain  gage  data  in  combination  with  finite  element  analyses.  The  stress  levels 
for  Zones  V,  VII-IX  were  determined  using  a coarse  grid  finite  element  analysis  and  a 
theoretical  stress  distribution  for  a circular  hole  in  an  infinite  plate  under  uniaxial 
tension.  The  stress  levels  for  Zones  VI  and  X were  determined  from  a fine  grid  finite 
clement  analysis. 

Fractographic  results  for  three  data  sets  (i.e.,  AFXLR4 , AFXMR4 , and  AFXHR4 ) for 
maximum  stress  levels  of  220.7  MPa  (32  ksi) , 234.5  MPa  (34  ksi) , and  262.1  MPa  (38  ksi) 
were  used  tc  calibrate  the  IFQ  model  parameters.  A 400-hour  block  load  spectrum  was 
used.  The  AFX  series  specimens  were  designed  for  15%  load  transfer.  The  specimens 
were  made  of  747S-T7351  aluminum  and  contained  two  MS90353-08  (h  dia.)  blind,  counter- 
sunk rivets  as  shown  m Fig.  7.  All  specimens  reflect  typical  aircraft  production 
quality,  tolerances  and  fastener  fits.  Nine  specimens  were  tested  per  data  set. 


STRESS 

ZONE 

LIMIT  STRESS  LEVEL,  MPa  (ksi) 

MJMBER  OF 
FASTENER  HOLES 

I 

195.2  (28.3) 

59 

II 

186.2  (27.0) 

320 

III 

167.6  (24.3) 

680 

IV 

115.2  (16.7) 

469 

V 

195.9  (28.4) 

8 

VI 

201.4  (-29.2) 

30 

V]  1 

223.5  (32.4) 

8 

VIII 

180.7  (26.2) 

8 

IX 

180.7  (26.2) 

12 

X 

177.2  (25.7) 

20 

1614 


330.2  mm  (13.0  ) 


meS 


S3I 


Figure  7 - IFQ  specimen 


The  crack  growth  rate  parameter  Q.  in  Eq.  2 (b  = 1)  was  determined  for  each  of  the 
three  AFX  data  sets.  Qq  was  determined  from  the  fractographic  results  using  a least- 
square  fit  of  Eq.  2 (Ret.  22).  A fractographic  crack  size  range  of  0 i27  mm  - 2.54  mm 
(0.005  inch  - 0.10  inch)  was  used.  An  upper  bound  EXFS  of  X = 0.762  mm  (0.03  inch) 
was  assumed  for  the  IFQ  distribution.  Using  Eq.  11  and  the  Estimated  Q„  values,  the 
corresponding  lower  bound  of  TTCI  value,  e,  for  each  reference  crack  size,  a0,  was 
determined  for  each  data  set.  The  results  of  Qq  and  e are  shown  in  Table  2. 

Table  2 - Summary  of  IFQ  model  parameters  for  fighter  spectrum 


DATA 

SET 

0"  MAX 
(MPa ) 

ao 

(mm) 

Qnxi04 

(HRS’1) 

e 

(MRS) 

ao 

60 

(HRS) 

0.76 

0 

15,033 

AFXLR4 

220.7 

1.27 

1.201 

4,253 

12,916 

2.54 

10,025 

13,421 

0.76 

0 

1.823 

8,721 

AFXMR4 

234. 5 

1.27 

2.037 

2,508 

7,759 

2. 54 

5,910 

9,093 

0.76 

0 

5,469 

AFXHR4 

262.1 

1.27 

4.731 

1,079 

5,098 

2.54 

2,545 

4 , r98 

1.777 

1.581 

1.852 


2.58'1 

2.412 

2.175 


POOLED  .-tQ  = 1.823,  AVERAGE  QQB0  = 1.928 


Fractographic  results  for  the  three  AFX  series  data  sets  were  combined  toaether 
to  determine  the  corresponding  "pooled"  a.  value  using  the  following  procedures. 
Time-to-crack-initiation  (TTCI)  results  for  three  different  reference  crack  sizes  (a„  = 
0.762  mm,  1.27  mm  and  2.54  mm)  were  used  for  each  of  the  three  data  sets.  The  adjusted 
ri’CI  data,  i.e.  , TTCI-  " data  for  each  reference  crack  size  for  each  data  set  were  normalized 
using  the  corresponding  average  values  (X) . Results  for  the  three  data  sets  were 
pooled  together  and  the  (TTCI-£  ) ,'X  data  were  ranked  in  ascending  order.  Equation  1 
was  transformed  into  a least-squares  fit  form  for  determining  the  pooled  un  v< lue  (Ref. 
22).  The  pooled  value  was  found  to  be  1.823  (Table  2).  0 

After  determining  «q  for  the  pooled  data  sets,  the  adjusted  TTCI's  for  each 
reference  crack  size  for  each  data  set  were  considered  separately  to  determine  the 
corresponding  6 values  (Ref.  22).  These  values  are  presented  in  Table  2.  Also 
summarized  m Table  2 are  the  Qq8q  values  for  the  nine  cases  considered.  For  generic 


I : 


EIFS , the  and  Q„  BQ  values  should  be  constants.  Average  values  of  Qq  8q  = 1.928  and 
Q - 1.823  are  usedJfor  the  present  durability  analysis.  A plot  of  Q versus  6 is 
shown  in  Fig.  3 for  the  three  data  sets  considered  (9  cases) . ' L' 


Figure  8 - Qq  versus  BQ  for  fighter  400-hour  load  spectrum 


The  crack  growth  rata  parameter  Q.  for  the  three  data  sets  is  plotted  against  the 
applicable  gross  stress  for  each  data  set  in  Fig.  9.  Using  a least-square  fit  (solid 
line  in  Figure  9) , the  following  expression  is  obtained  for  Q as  a function  of  stress 
level  when  stress  is  expressed  m ksi  units: 

Qi  = 1.427  X 10~16  a7*928  (20) 

When  stress  is  expressed  m MPa  units,  the  appropriate  expression  for  fk  is  as  follows: 


3.2  X 10 


-23  7.928 


(21) 


Equation  20  is  used  to  estimate  the  value  for  each  of  the  ten  stress  regions 
shown  in  Fig.  6. 

Crack  exceedance  predictions  for  the  fighter  lower  wing  skin  were  determined  using 
Eqs..  13,  15,  and  20  as  well  as  the  following  parameters:  x = 0.762  m (0.03  inch), 

a = 1.823,  Q 3 = 1.928  (average),  c = 4 and  various  t values.  The  results  are  pres- 

ented in  various  formats  as  described  below. 


The  extent  of  damage  predictions  for  the  fighter  lower  wing  skin  are  summarized  in 
Table  3 at  t = 16000  flight  hours  for  each  of  the  ten  stress  regions  shown  in  Fig.  6. 

The  number  of  fastener  holes  with  a crack  si/e  >_0.762  mm  (0.03  inch),  L(x),  and  the 
standard  deviation , uT  ( i ) , was  estimated  to  be  17.6  and  4.077,  respectively.  Based  on 
the  test  results  for"the  right  hand  and  left  hand  lower  wing  skins,  an  average  of  16.5 
fastener  holes  had  a crack  size  ^0.762  mm  (0.03  inch)  all-  16000  hours.  In  Table  3, 
the  predicted  extent  of  damage  results  track  the  average  test  results  for  the  individual 
stress  regions  very  well. 

In  Fig.  10  the  predicted  percentages  of  crack  exceedance  versus  fastener  hole  crack 
size  are  plotted  for  the  fighter  lower  wing  skin  at  t = 16000  flight  hours.  Curves  1,  2 
and  3 are  based_on  L(t)  x 100%/N*,  (L(t)  + aT(x)]  x 100%/N*  and  (L(t)  - a (x)]  x 100%/N*, 
respectively.  L(x)  and  oT  (x)  are  defined  by  tqs.  18  and  19,  respectively.  N*  is  the 
total  number  of  fastener  holes  m the  fighter  lower  wing  skin  (i.e.,  1614  holes).  Since 
the  number  of  fastener  holes  m each  stress  region  is  large,  it  is  reasonable  to  approx- 
imate the  binomial  distribution  by  the  normal  distriVmtion.  The  corresponding  exceedance 
probabilities  for  curves  1,  2,  3 are  shown  in  Fig.  10  m parentheses. 
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Figure  10  - Percentage  of  crack  exceedance  versus  crack  size  at  16,000 
hours  for  3 probability  levels  (fighter) 

Test  results  tor  the  right  and  left  hand  lower  wing  skin  (at  X.  = 0.762  mm  and  t = 
16000  hoUi.s)  are  plotted  as  a circle  and  a square,  respectively,  inrig.  10.  Approxi- 
mately 1.1%  of  the  fastener  holes  m the  fighter  lower  wing  skin  are  predicted  to  have 
a crack  size  >_  0.762  mm  (0.03  inch)  at  t - 16000  hours..  This  compares  'with  an  average 
of  1.02%  based  on  test  results  for  the  right  hand  and  left  hand  lower  wing  skins. 

In  Fig.  10,  the  predicted  average  percentage  of  crack  exceedance  decreases  rapidly 
for  larger  crack  sizes.  For  example,  the  average  percentage  of  crack  exceedance  for 
the  fighter  lower  wing  skin  decreases  from  approximately  1.1%  at  = 0.762  mm  (0.03  ") 
to  approximately  0.14%  at  Xj  = 1.27  mm  (0.05  inch).  Crack  exceedance  predictions  are 
based  on  the  service  crack  growth  master  curve  defined  by  Eqs.  13  and  14.  A single 
service  crack  growth  master  curve  may  not  adequately  fit  the  full  range  of  desired 
crack  sizes  for  all  crack  exceedance  predictions.  For  example,  different  service  crack 
growth  master  curves  are  required  to  fit  two  different  crack  size  ranges  as  illustrated 
m Fig.  11.  Curve  1 and  Curve  2 shown  in  Fig.  11  apply  to  crack  size  ranges  Aj  and  A~, 
respectively.  Crack  exceedance  predictions  based  on  Curves  1 and  2 of  Fig.  11  will  be 
different  for  the  same  crack  exceedance  size,  X . For  example,  p(i,T)  predictions 
based  on  Curve  2 for  X^,  and  Xj,  t2  will  bexlarger  than  those  based  on  Curve  1. 

The  extrapolation  of  crack  exceedance  predictions  to  larger  crack  sizes  s houid  be 
consistent  with  the  applicable  ciack  growth  process  for  given  design  conditions  und  the 
crack  exceedance  crack  size,  X^.  Further  research  is  needed  to  develop  a better 
understanding  and  confidence  m crack  exceedance  predictions  for  different  crack  sizes, 
materials,  and  design  conditions. 


Analytical  predictions  of  the  extent  of  damage  are  presented  in  Fig.  12  in  an 
exceedance  probability  format.  In  this  case,  the  predicted  number  of  fastener  holes  in 
the  fighter  lower  wing  skin  with  a crack  size  > 0.762  mm  (0.03  inch)  are  plotted  as  a 
function  of  flight  hours  for  different  exceedance  probability  values  (i.e.,  P = 0.05, 
0.50,  0.95).  The  plots  are  based  on  Eq.  15,  x,  = 0.762  mm  (0.03  inch),  i.  = 1.823,  Q„B 

= 1.928  (average),  f=  4,  N = 1 614  fastener  holes,  Z - ii.65  and  L(x)  iZo^(t1.  For 

example,  at  " - 16,000  houJs,  L ( t ) = 17.6  fastener  holes  and  ut(t)  = 4.07<.  The  upper 
bound  prediction, L ( i + Zu.ii),  is  approximately  24.3  fastcner~holes . In  other  words, 
there  is  a probability  of  Cf . 0 5 that  more  than  24.3  fastener  holes  m the  fighter  lower 
wing  skin  will  have  a crack  size  > 0.762  mm  (0.03  inch)  at  16000  flight  hours.  There 
is  a probability  of  0.50  and  0.95,  respectively,  that  more  than  17.6  and  10.9  fastener 
holes  will  have  a crack  size  >0.762  nun  (0.03  inch)  at  ' = 16000  flight  hours.  The 


NUMBER  OF  HOLES  WITH  CRACK 


stress  zone  were  all  increased  by  the  same  percentage.  The  results  shown  m Fig.  13 
can  be  used  to  assess  the  extent  of  damage  as  a function  of  stress  level  and  flight 
hours.  This  format  is  particularly  useful  for  evaluating  durability  design  tradeoffs 
m terms  of  the  extent  of  damage.  For  example,  at  t = 16000  flight  hours  approximately 
1.1%  of  the  fastener  holes  in  the  fighter  lower  wing  skin  would  be  predicted  to  exceed 
a crack  size  of  0.762  mm  (0.03  inch)  for  the  baseline  stress  levels.  If  the  baseline 
stresses  were  increased  to  1.1a  and  i.2a,  the  predicted  average  percentage  of  holes 
with  a crack  size  >_  0.762  mm  (0.03  inch)  would  be  approximately  4%  and  12%,  respectively. 
This  provides  a quantitative  measure  of  the  structural  durability  as  a function  of 
stress  level  and  flight  hours. 


SYMBOLS 


Figure  13  - Average  percentage  of  holes  with  crack  size  > 0.76  mm 
(0.03  inch)  versus  flight  hours  - stress  level  format 
( f lahter ) 


fi.  Comp iex-bpi ice  Specimens  Subjected  to  Bomber  Load  Spectrum 

A durability  analysis  of  complex-splice  specimens  subjected  to  a bomber  load  spectrum 
is  presented.  *»nalytical  predictions  of  the  except  of  damage  m the  specimens  are 
presented  m variOub  formats  and  compared  with  f ractographic  resuits.  The  analytical/ 
experimental  results  are  summarized  here  and  described  m more  detail  m Ref.  27. 


The  complex-splice  specimen  geometry  is  presented  m Fig.  14.  The  specimens  were 
made  ci  7475-T735I  aluminum  plate  and  countersunk  steel  rivots.  A bomber  load  spectrum 
was  applied.  Based  on  a simplified  stress  analysis  and  strain  gaae  results,  the 
maximum  gross  stroas  ip.  the  outer  row  of  fastener  holes  at  the  faying  surface  was 


MATERIAL:  7475-T7351  ALUMINUM 


Figure  14  - Complex  --  splice  specimen 


After  testing,  ail  fastener  holes  m the  outer  rows  were  inspected.  Fractography 
was  performed  for  the  largest  crack  in  each  fastener  hole  in  the  outer  rows.  Twenty 
five  out  of  110  fastener  holes  m the  outer  rows  had  a crack  size  >_1 . 27  mm  (.05  inch) 
at  13,500  hours,  Hence,  1:2.7%  of  the  fastener  holes  in  the  outer  rows  had  a crack  size 
H.27  mm  (.05  inch)  at  13,500  hours. 


The  IF,_>  of  the  fastener  holes  was  based  on  the  fractngraphic  results  for  nine  data 
sets  and  three  different  reference  crack  sizes.  The  specimens  were  made  of  7475-T7351 
aluminum  and  contained  ? countersunk  rivets,  toad  transfe*  levels  of  i5%,  30%  ..no  40% 
were  considered.  All  specimens  had  the  same  configuration  dig.  7}  with  the  same 
overall  length  and  basic  test  section  dimensions.  However,  the  lug  end  dimensions 
varied  depending  on  the  amount  of  load  transfer.  Three  maximum  stress  levels  were 
considered  for  each  load  transfer  leve>. 


A f ractographic  crack  size  range  of  C.127  mm  - 2.54  ran,  (0.005  inch  - 0.1  inch)  was 
considered,  An  upper  bound  dIFS  of  X = 1,27  r«m  (0.05  inch)  was  assumed  for  the  IFQ 
distribution.  A fractcgraphy  scalmg"factor  oi  i = 4 was  used.  Tne  same  data  pooling 
procedures  wo i e used  winch  were  previously  described  for  the  lighter  desionstration . 

The  average  a 


and  Qnpn  values  were  found  to  be  2.702  and  2.823,  respectively, 


£> 


m 

Nsa 


/ 0 - / *i>  The  crack  growth  rate  parameter  for  the  9 data  sets  is  plotted  against 

the  applicable  gross  stress  for  each  data  set  in  Fig.  15.  The  solid  line  represents 
the  least-square  best  fit  through  the  plot  points.  The  dashed  lines  have  the  same 
slope  as  the  solid  line  and  they  encompass  all  the  plot  points.  The  corresponding 
best-fit  equation  for  as  a function  of  gross  stress  level  when  stress  is  expressed 
in  ksi  units  is  as  follows: 

Q.  = 6.151  X Itf13  a5’331  <22> 

When  stress  is  expressed  in  MPa  units,  the  appropriate  expression  for  Q is  as  follows: 


Qi  = 1.895  X 10~17  o5'381 
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Qt  = 1.895X10'1  JrS.381 
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Figure  15  - Q.  versus  gross  stress  for  bomber  load  spectrum 


Crack  exceedance  predictions  for  the  complex-splice  specimens  were  determined 
using  Eqs.  13,  15,  and  22.  At  x = 13,500  hours,  an  average  of  9 fastener  holes  (8.3%) 
were  predicted  to  exceed  a crack  size  of  1.27  nun  (0.05  inch).  The  test  results  showed 
an  average  or  25  fastener  holes  (22.7%)  exceeding  a crack  size  of  1.27  mm  (0.05  inch). 
The  difference  in  the  predicted  and  test  clock  exceedances  is  attributed  mainly  to  the 
stress  level  used  m the  predictions.  The  actual  stress  level  and  distribution  m the 
outer  row  of  fastener  holes  is  far  more  complex,  due  to  lateral  bending  effects,  than 
those  considered  for  the  damage  assessment.  The  crack  exceedance  predictions  are  very 
sensitive  to  the  gross  applied  stress  level  used.  This  is  illustrated  in  Fig.  16.  The 
solid  line  represents  average  crack  exceedance  predictions  for  the  gross  stress  level 
of  246.8  MPa  (35.8  ksi)  obtained  using  the  simplified  stress  analysis  approach.  The 
dashed  lines  repiesent  average  crack  exceedance  predictions  for  other  gross  stress 
levels.  Also  plotted  as  a single  point  is  the  average  test  crack  exceedance  at  i = 
13,500  hours.  It  can  be  seen  that  if  the  gross  applied  stress  level  used  in  the 
predictions  were  266.1  MPa  (38.6  ksi)  rather  than  246.8  MPa  (35.8  ksi),  the  predicted 
crack  exceedance  at  x=  13,500  hours  would  match  the  test  results.  Hence,  a more 
accurate  stress  analysis  could  result  m improved  predictions. 


XnAiOGi.’i *-2  C*. 
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Figure  16  - Average  percentage  holes  with  crack  size  > 1.27  mm 

(0.05  inch)  versus  flight  hours  stress  level  format 
( bomber ) 

Other  useful  crack  exceedance  formats,  previously  discussed  for  the  fighter  demon- 
stration, are  presented  in  Figs.  17  and  18  for  tne  complex-splice  specimens. 

VI  CONCLUSIONS  AND  DISCUSSIONS 


Piobabilistic  fracture  mechanics  methods  for  durability  analysis  have  been  described 
and  demonstrated  for  both  a full-scale  fighter  aircraft  structure  and  for  a complex 
splice  subjected  to  a bomber  spectrum.  These  methods  can  be  used  to  analytically  assure 
compliance  with  the  Air  Force's  durability  design  requirements.  The  analytical  tools 
described  can  be  used  to  quantify  the  extent  of  damage  as  a function  of  the  durability 
design  variables  for  structural  details  m a part,  a component  or  airframe.  Once  the 
economic  life  and  durability  critical  parts  criteria  are  established,  the  extent  of 
damage  predictions  can  be  used  to  assure  design  compliance  with  Air  Force  durability 
requirements . 

An  initial  fatigue  quality  model  can  be  used  to  define  the  EIFS  cumulative  distri- 
bution using  suitable  fractographic  results.  Procedures  and  guidelines  have  been 
developed  for  determining  the  IFQ  model  parameters  for  pooled  fractographic  data  sets 
and  for  scaling  TTCI  results.  The  parameters  and  QqSq  provide  the  basis  for  putting 
fractographic  results  on  a common  baseline  for  quantifying  the  initial  fatigue  quality. 
For  generic  EIFS,  a.  and  QgE-  should  be  constants  for  different  fractographic  data  sets 
(same  material,  fastener  type/fit,  and  drilling  technique),  loading  spectra,  stress 
levels  and  percent  load  transfer.  Encouraging  results  have  been  obtained  to  justify 
the  use  of  the  same  EIFS  cumulative  distribution  for  crack  exceedance  predictions  for 
different  design  conditions.  Further  research  is  required  to  confirm  the  IFQ  distri- 
butions for  different  materials,  load  spectra,  stress  levels,  fastener  types/diameters/ 
fit,  % load  transfer,  etc.  A considerable  amount  of  fractographic  results  exist  which 
need  to  be  evaluated  using  the  IFQ  model. 

The  effects  of  fretting,  clamp-up,  corrosion,  size  effect  (scale-up  from  coupon  to 
component),  faying  surface  sealant,  interference-fit  fasteners,  etc.  on  IFQ  need  to  be 
investigated.  Also,  the  feasibility  of  using  no-load  transfer  specimens  with  multiple 
holes  for  quantifying  the  IFQ  should  be  evaluated  using  spectrum  and  constant  amplitude 
loading.  This  could  provide  an  economical  way  to  generate  the  fractographic  results 
needed  to  quantify  the  IFQ. 


Theoretically,  the  IFQ  model  can  be  used  to  quantify  the  EIFS  cumulative  distri- 
bution for  various  structural  details  as  long  as  fractographic  results  are  available 
for  the  details  to  be  included  in  the  durability  analysis.  The  IFQ  model  has  been 
evaluated  using  fractographic  results  for  fastener  holes.  Suitable  specimens  and 
guidelines  need  to  be  developed  for  generating  crack  initiation  and  crack  growtn 
results  for  other  details  such  as,  cutouts,  fillets,  lugs,  etc.  Fractographic  results 
should  be  developed  and  evaluated  for  such  details  so  that  tne  durability  analysis 
methods  described  can  be  efficiently  applied  to  different  types  of  structural  details 
in  typical  aircraft  structures. 

The  accuracy  of  crack  exceedance  predictions,  based  on  the  same  EIFS  cumulative 
distribution,  needs  to  be  evaluated  for  different  design  conditions.  Also,  IFQ  model 
parameter  sensitivity  studies  need  to  be  performed  to  better  understand  the  average 
parameter  values  and  variances  and  the  impact  of  these  parameters  on  the  IFQ  for 
different  fractographic  data  sets. 

The  durability  analysis  methodology  was  developed  for  crack  exceedance  predictions 
for  relatively  small  crack  sizes  (e.g.,  = 2.54  nan)  m structural  details.  The  largest 
crack  in  each  detail  was  assumed  to  be  statistically  independent  to  justify  using  the 
binomial  distribution  for  combining  crack  exceedance  predictions  for  structural  details. 
If  the  largest  crack  m a given  detail  doesn't  significantly  affect  the  growth  of 
cracks  in  neighboring  details,  perhaps  the  proposed  durability  analysis  methodology  can 
be  extended  to  crack  s^zes  >2.54  mm  (0.10  inch).  The  simplistic  crack  growth  rate 
equation  (Eq.  12)  is  not  suitable  for  use  m the  crack  exceedance  predictions  for  crack 
sizes>2.54  mm  (0.30  inch) . However,  a general  service  crack  growth  master  curve  can  be 
generated  under  given  design  conditions  which  is  valid  for  crack  sizes  >2.54  mm  (Ref. 

11,  12  15,  and  22).  Nevertheless,  this  approach  has  not  been  demonstrated  m the 
present  study  and  further  research  is  required  to  extend  the  probabilistic,  fracture 
mechanics  approach  developed  to  larger  crack  sizes. 

Two  different  Fa/0)(x)  equations  (i.e.,  Eqs.  7 and  9)  were  presented  for  represent- 
ing the  IFQ.  Eitherequation  works  but  Eq.  9 is  recommended  for  two  reasons:  (1)  It 

assures  all  EIFS'sin  the  IFQ  distribution  will  be>0,  and  (2)  the  crack  growth  rate 
parameter  Q0  can  be  easily  determined  from  the  fractographic  results  and  the  resulting 
Qn  values  for  different  data  sets  will  be  directly  comparable.  If  Eq.  7 is  used,  a 
common  b parameter  (Eq,  2)  must  be  imposed  for  different  fractographic  data  sets  to  put 
the  Q.fi.  values  on  a comparable  baseline.  As  long  as  b > 1,  all  EIFS's  in  Eq.  7 will 
be  > o.  Further  studies  are  needed  to  evaluate  the  accuracy  of  these  two  F ...  (x) 

equations. 


The  EIFS  cumulative  distribution,  F ... (x),  is  independent  of  the  reference  crack 
size,  a..  This  is  illustrated  in  Eqs.  7?an<3  9.  Therefore,  the  TTCI  distribution  for 

different  reference  crack  sizes  will  transform  into  a common  F ... (x) . 

a luj 

The  IFQ  model  is  simply  a "mathematical  tool"  for  quantifying  the  IFQ  of  structural 
details.  Therefore,  the  resulting  EIFS's  must  be  considered  in  the  context  of  the  IFQ 
model  and  the  fractographic  results  used  to  calibrate  the  model  parameters.  EIFS's 
should  be  considered  as  hypothetical  cracks  used  for  crack  exceedance  predictions 
rather  than  actual  initial  flaws  per  se. 

Back  extrapolations  of  fractographic  data  must  be  done  consistently  to  put  the 
EIFS's  on  a common  baseline  for  different  data  sets.  Inconsistent  EIFS  results  will  be 
obtained  if  the  EIFS  distribution  is  determined  by  back  extrapolating  the  fractcgraphy 
results  for  individual  specimens  and  then  fitting  a statistical  distribution  to  the 
EIFS  results  for  different  data  sets.  Two  problems  result  n this  approach  is  used: 

(1)  the  EIFS's  are  not  on  a common  baseline  for  different  data  sets,  and  (2)  the 
result Lng  EIFS  distribution  is  not  statistically  compatible  t ith  the  TTCI  distribution 
and  the  fatigue  wear  out  process.  The  resulting  EIFS  distribution  should  be 
statistically  compatible  with  the  TTCI  distribution.  The  IFQ  model  presented  in  this 
paper  satisfies  this  requirement. 

Several  useful  applications  of  the  durability  analysis  methodology  developed  are: 

(1)  The  evaluation  of  durability  design  tradeoffs  in  terms  of  structural  design  variables, 

(2)  the  evaluation  of  structural  maintenance  requirements  before  or  after  aircraft  is 
committed  to  service,  and  (3)  the  evaluation  of  aircraft  user  options  affecting  life-cycle- 
costs,  structural  maintenance  requirements,  and  operational  readiness. 
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